Phantom dark energy with varying-mass dark matter particles: acceleration and 

cosmic coincidence problem 
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We investigate several varying-mass dark-matter particle models in the framework of phantom 
cosmology. We examine whether there exist late-time cosmological solutions, corresponding to 
an accelerating universe and possessing dark energy and dark matter densities of the same order. 
Imposing exponential or power-law potentials and exponential or power-law mass dependence, we 
conclude that the coincidence problem cannot be solved or even alleviated. Thus, if dark energy 
is attributed to the phantom paradigm, varying-mass dark matter models cannot fulfill the basic 
requirement that led to their construction. 
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I. INTRODUCTION 

Recent cosmological observations support that the uni- 
verse is experiencing an accelerated expansion, and that 
the transition to the accelerated phase realized in the re- 
cent cosmological past [l| . In order to explain this unex- 
pected behavior, one can modify the theory of gravity 
or introduce the concept of dark energy which provides 
the acceleration mechanism. The most explored dynami- 
cal dark energy models of the literature consider a canon- 
ical scalar field (quintessence) a phantom field, that 
is a scalar field with a negative sign of the kinetic term 
or the combination of quintessence and phantom in 
a unified model named quintom Q. 

The dynamical nature of dark energy introduces a new 
cosmological problem, namely why are the densities of 
vacuum energy and dark matter nearly equal today al- 
though they scale independently during the expansion 
history. The elaboration of this "coincidence" problem 
led to the consideration of generalized versions of the 
aforementioned scenarios with the inclusion of a coupling 
between dark energy and dark matter. Thus, various 
forms of "interacting" dark energy models [1, 0, B Q 
have been constructed in order to fulfil the observational 
requirements. In the case of interacting quintessence one 
can find accelerated attractors which moreover give dark 
matter and dark energy density parameters of the same 
order, thus solving the coincidence problem but 
paying the price of introducing new problems such is the 
justification of a non-trivial, almost tuned, sequence of 
cosmological epochs . In interacting phantom models 
B B [3 J ths existing literature remains in some spe- 
cial coupling forms which suggest that the coincidence 
problem might be alleviated 

An equivalent approach is to assume that dark energy 
and dark matter sectors interact in such a way that the 



dark matter particles acquire a varying mass, dependent 
on the scalar field which reproduces dark energy p^ . 
This consideration allows for a better theoretical justi- 
fication, since a scalar-field-dependent varying-mass can 
arise from string or scalar- tensor theories [T^. Indeed, 
in such higher dimensional frameworks one can formu- 
late both the appearance of the scalar field (which is 
related to the dilaton and moduli fields) and its effect 
on matter particle masses (determined by string dynam- 
ics, supersymmetry breaking, and the compactification 
mechanism) [l6j . In quintessence scenario, such varying- 
mass dark matter models have been explored in cases 
of linearM, Hi, [H, [13, power-law [l^ or exponential 
0,113, nil scalar- field dependence. The exponential case 
is the most interesting since, apart from solving the co- 
incidence problem, it allows for stable scaling behavior, 
that is for a large class of initial conditions the cosmo- 
logical evolution converges to a common solution at late 
times M,M. 
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In the present work we are interested in investigating 
varying-mass dark matter models in scenarios where dark 
energy is attributed to a phantom field. Although such 
a framework could lead to instabilities at the quantum 
level [2^ , there have been serious attempts in overcoming 
these difficulties and construct a phantom theory consis- 
tent with the basic requirements of quantum field theory, 
with the phantom fields arising as an effective description 
[23| . Performing a complete phase-space analysis using 
various forms of mass-dependence and scalar-field poten- 
tials, we examine whether there exist stable late-time ac- 
celerating solutions which moreover solve the coincidence 
problem. As we will show, the coincidence problem can- 
not be solved in any of the investigated models. 

The plan of the work is as follows: In section [ll] we 
construct varying-mass dark matter models in the frame- 
work of phantom cosmological scenario and we present 
the formalism for the transformation into an autonomous 
dynamical system. In section Hill we perform the phase- 
space stability analysis for four different models, using 
various mass-dependence forms and phantom potentials, 
and in section IIVI we discuss the corresponding cosmo- 
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logical implications. Finally, in section |V] we summarize 
the obtained results. 



II. VARYING-MASS DARK MATTER 
PARTICLES IN THE FRAMEWORK OF 
PHANTOM COSMOLOGY 

Let us construct a cosmological model where dark en- 
ergy is attributed to a phantom field, in which the dark 
matter particles have a varying mass depending on this 
field. Throughout the work we consider a flat Robertson- 
Walker metric: 



dt^ ~ a^(t)dy? 



(1) 



with a the scale factor and t the comoving time. 

In the phantom cosmological paradigm the energy den- 
sity and pressure of the phantom scalar field (j) are: 



P0 



0' + 1^(0) 



(2) 
(3) 



Since general covariance leads to total energy conser- 
vation, we deduce that the evolution equation for the 
phantom energy density will be: 



3H{p^+p^) = - 



1 



dMDM{<l>) 



MdkM) d(j) 



4>Pdm- (8) 



Thus, 



< corresponds to energy transfer from 



dark matter to dark energy, while '^^^'^'^'"'^^ 4> > corre- 
sponds to dark energy transformation into dark matter. 

Equivalently, using the definitions ([2]) and the 
phantom evolution equation can be written in field terms 
as: 



. -t- 3H<j> - 



dV{(f>) 



1 



dMDM{4') 



■ RDM- 



(9) 



Finally, the system of equations closes by considering the 
Fricdmann equations: 



-TiP't> + Pdm), 



(10) 



where V{(f>) is the phantom potential and the dot denotes 
differentiation with respect to comoving time. In such a 
scenario, the dark energy is attributed to the phantom 
field, and its equation of state is given by 



WDE 



(4) 



As was mentioned in the introduction, in varying-mass 
dark matter models the central assumption is that the 
dark-matter particles have a (/)-dependent mass M^m 
while dark matter is considered as dust. Thus, for the 
dark matter energy density we have the standard defini- 
tion 



PDM = MDMi4')nDM, 



(5) 



where is the number density of the dark-matter 

particles. As usual, in the case of FRW geometry, it is 
determined by the equation 



riDM + 'SHuDM = 0, 



(6) 



with H the Hubble parameter. Therefore, differentiating 
^ and using ^ we obtain the evolution equation for 
PDM, namely: 



Pdm + 'iHpDM = -T-) 



1 



dMDM{<P) 



Mdm{4') d(j) 



4>Pdm- (7) 



Obviously, in a case of (/)-independent dark-matter par- 
ticle mass, we re-obtain the usual evolution equation 
Pdm + iHpDM = 0. Therefore, we observe that the 
0-dependent mass reveals the interaction between dark 
matter and dark energy (that is the phantom field) sec- 
tors that lies behind it. 



H = — —[P4,+P4, + PDM 



(11) 



where we have set = SttG. Although we could 
straightforwardly include baryonic matter and radiation 
in the model, for simplicity reasons we neglect them. 

Alternatively, one could construct the equivalent un- 
coupled model described by: 



Pdm 



where 



WDM,eff 



3H{1 + WDM,eff)pDM = 
9^ + 3H{1 + W^^eff)P4> = 0, 



dMDM{4>) _±_ 

3H 



MDM{(t>) dc^ 

1 



(12) 
(13) 



(14) 



dMDMicj}) (j) Pdm 

^ • (15) 



Mdm{(I)) d(j) 3H pcjy 



However, it is more convenient to introduce the "total" 
energy density ptot = Pdm + P<f,, obtaining: 



ptot + 3H{1 + Wtot)ptot = 0, 



with 



Wtot 



p<ji + Pdm 



(16) 



(17) 



where ft^ = = floE- Obviously, since ptot — 

3H^/k'^, leads to a scale factor evolution of the form 
a(t) (X t^/i^i^+^tot)) ^ jj^ ^jjg constant Wtot case. However, 
at the late-time stationary solutions that we are study- 
ing in the present work, wtot has reached to a constant 
value and thus the above behavior is valid. Therefore, we 
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conclude that in such stationary solutions the condition 
for acceleration is just Wtot < —1/3. 

In order to perform the phase-space and stability anal- 
ysis of the phantom model at hand, we have to trans- 
form the aforementioned dynamical system into its au- 
tonomous form [23, HBl • This will be achieved by intro- 
ducing the auxiliary variables: 



K(j) 



y = 



K(j) 



(18) 



together with AI — In a. Thus, it is easy to see that for 



every quantity F we acquire F 
variables we obtain: 



H 



dF 
dM- 



„ _ K P0 2,2 



Using these 



(19) 



and 



W4, 



Wtot 



y 



(20) 



(21) 



We mention that relations ([^0]) and pT|) are always valid, 
that is independently of the specific state of the system 
(they are valid in the whole phase-space and not only 
at the critical points). Finally, note that in the case of 
complete dark energy domination, that is pdm and 
ri^ — > 1, we acquire Wtot ~ Wcf, < —1, as expected to 
happen in phantom-dominated cosmology. 

The next step is the introduction of a specific ansatz 
for the phantom potential V(0), and a specific ansatz for 
the dark- matter particle mass function M£)m(0)- In this 
case the equations of motion (UHl) and pT|) can 

be transformed into an autonomous system containing 
the variables x and y and perhaps z {z is present only for 
some such ansatzes) and their derivatives with respect to 
M = Ina. 

Having transformed the cosmological system into its 
autonomous form: 



X' = f(X), 



(22) 



where X is the column vector constituted by the auxil- 
iary variables, f(X) the corresponding column vector of 
the autonomous equations, and prime denotes derivative 
with respect to M = Ina, we extract its critical points 
Xc satisfying X' = 0. Then, in order to determine the 
stability properties of these critical points, we expand 
(22) around Xc , setting X = Xc -I- U with U the pertur- 
bations of the variables considered as a column vector. 



Thus, for each critical point we expand the equations for 
the perturbations up to the first order as: 



U' = Q U, 



(23) 



where the matrix Q contains the coefficients of the per- 
turbation equations. Thus, for each critical point, the 
eigenvalues of Q determine its type and stability. 



III. PHASE-SPACE ANALYSIS 

In the previous section we constructed a cosmological 
scenario where the dark matter particles have a varying 
mass, depending on the phantom field. Additionally, we 
presented the formalism for its transformation into an 
autonomous dynamical system, suitable for a stability 
analysis. In this section we introduce specific forms for 
V{(l)) and MoMi4')j and we perform a complete phase- 
space analysis. 

For the scalar field potential we consider two well 
studied cases of the literature, namely the exponential 
MM: 



V{(j)) = yoe""^'"^ 
and the power-law one [H, [2^ : 



(24) 



(25) 



For the dark matter particle mass we consider two possi- 
ble cases, namely an exponential dependence [l9l [20l [2]| : 

AfoM(0) = Moe-^^'^ (26) 

and the power-law one [l8| : 

Mdm{<I)) Mo^-^\ (27) 

Therefore, in the following we consider four different 
models, arising from the aforementioned combinations. 



A. Model 1: Exponential potential and 
exponentially-dependent dark-matter particle mass 



Inserting the auxiliary variables (|18p into the equations 
of motion Q, (H]), ^U\j and ^T^j, we result in the follow- 
ing autonomous system: 



a; = -3a; + ^a;(l - - y^) - J 



|Aii(l + 2;'-y') 



3 3 

y' = ^'^{^ -x^ -y"^)- \j 2^^! ^y- 



(28) 



Note that in this case, the auxiliary variable z is not 
needed. 
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The critical points {xc,yc) of the autonomous system 
are obtained by setting the left hand sides of the 
equations to zero. The real and physically meaningful 
(that is corresponding to y > and < il^ < 1) of them 
are: 



Q 



Xc2 = 



Ai - /ii 



yc2 



Ml (Ai - Ml) 



|Ai - Mil 



(29) 



and in table |T] we present the necessary conditions for 
their existence. The 2x2 matrix Q of the linearized 
perturbation equations writes: 



i {-9x1 - - 3 (2/2 + 1)) 

-ij/c + \/6Ai) 



Vc (V6(mi - Ai) - 3xc) 
i {~9y^,-x, (3a;, + V6Ai) +3) 



Therefore, for each critical point of table IH we examine 
the signs of the real parts of the eigenvalues of Q, which 
determine the type and stability of this specific critical 
point. In table [T] we present the results of the stability 
analysis. In addition, for each critical point we calculate 
the values of Wfot (given by relation (|2ip ). and of 
(given by (|19p ). Thus, knowing wtot we can express the 
acceleration condition wtot < —1/3 in terms of the model 
parameters. 

The critical point A exists always and it is either a 
saddle point (the Q- eigenvalues have real parts of differ- 
ent sign) or an attractor (the Q-eigenvalues have nega- 
tive real parts). The critical point B, if it exists, it is 
always a saddle point. The cosmological model at hand 
admits another critical point, namely C, which is unphys- 
ical since it leads to fi^ < 0. This point has coordinates 

(^Xc3 = — -y/lMij Ucs = 0^ and it is either a saddle point 
or an attractor. If Mi(Mi — Ai) > 3/2 it is an attractor 
and in this case, although unphysical, it can attract an 
open set of orbits from the interior of the physical region 
of the phase space. 

In order to present this behavior more transparently, 
we evolve the autonomous system (1281) numerically for 
the parameters Ai — 0.4 and /ii = 2, and the results 
are shown in figure [1] Depending on which region of the 
phase-space does the system initiates, it lies in the basin 
of attraction of either A or C, and thus it is attracted 
by one or the other point. In particular, the orbits ini- 
tially below the stable manifold of B-points converge to- 
wards C, while the orbits initially above this curve con- 
verge to A. Interestingly, A is not the global attractor 
for points at the physical region (region corresponding 
to < 51^ < 1, bounded by the dashed (red) curves). 



However, if ^ 



2 



< < 



+ 



2 



point C is 



always a saddle one and B does not exist. Thus, in this 
case A is the attractor for all the points located at the 
physical region. This behavior is presented in figure [2l 
Finally, for completeness we mention that in the trivial 
case fii — the origin is also a saddle point. It repre- 
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sents matter-dominated universe {floM = sh^ ' ~ 1) 
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FIG. 1: (Color Online) Pftase plane of Model 1 for the pa- 
rameter values Ai — 0.4 and Hi — 2. The stable manifold of 
B (thick curve) divides the physical part of the phase space 
(regton corresponding to < Q,^ < 1, bounded by the dashed 
(red) curves) in two regions. The orbits initially below this 
curve converge towards C. The orbits initially above this curve 
converge to A. 



with ^-independent dark matter particle mass. 



B. Model 2: Power-law potential and 
power-law-dependent dark-matter particle mass 



Inserting the auxiliary variables ([T5)) into the equations 
of motion (l?]), dH), (dH]) and dTT]), we result in the follow- 
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Cr. P. 


Xc 


yc 


Existence 


Stable for 




Wtot 


Acceleration 


A 
B 


Xc2 


yc2 


Always 
min{/ii — 3, Ai + 3} > Ai/^i, 
Ml / -^1 


Ai (^1 - Ai) < 3 
Never 


1 


-1(3 + A?) 

Ml 


Always 
^ii < 0, /ii < Ai < —2/^1 
^ii > 0, — 2/ii < Ai < Ml 


(Ai-Mi)^ 


^1-Ml 



TABLE I: The real and physically meaningful critical points of Model 1 and their behavior. 




{Xcb = 0, ?/c5 = 1, 2c5 = 0) , (31) 

FIG. 2: (Color Online)P/iase plane of Model 1 for the param- 
eter values Ai = 1 and = 0.5. In this case the critical point 
B does not exist and all orbits initially at the physical region 

converge to A. The dashed (red) curves hound the physical and in table |TT] we present the necessary conditions for 
part of the phase space, that is corresponding toQ<Q,^<l. ^^^eir existence. The 3x3 matrix Q of the linearized 

perturbation equations writes: 



3) ycizcip2 - A2) - 3xe) i (2/2(^2 - A2) - (a;2 + 1) ^2) 

—2XrZr 



In this case, the critical points are non-hyperbohc, that 
is there exists always at least a zero eigenvalue. We men- 
tion that for non-hyperbolic critical points the result of 
linearization cannot be applied in order to investigate 
the local stability of the system (the system can be un- 
stable to small perturbations on the initial condition or 
to small perturbations on the parameters) [13, [H, [29j . 
However, it is possible to get information about the ex- 
istence and the dimensionality of the stable manifold by 
applying the center manifold theorem Doing so we 
deduce that the dimensionality of the local stable man- 
ifold is 1 and 2 for D and E respectively. In particular, 
the stable manifold of D is tangent, at the critical point, 
to the X-axis, while the stable manifold of E is tangent, 



at the critical point, to the xy-\Aa,Tie. The existence of an 
ID stable manifold for D, implies that the orbits asymp- 
totic to D as t ^ — 00 are contained in either an unstable 
or center manifold (each one of dimensionality 1, that is 
a curve). There are some exceptional orbits converging 
to D as i ^ -t-00, but these have a zero measure. On 
the other hand, the fact that E has a 2D stable manifold 
implies that there exists a non-zero-measure set of orbits 
that converges to E as t — > +00. Finally, there are some 
exceptional orbits contained in its center manifold that 
cannot be classified by linearization. In summary, using 
more sophisticated tools such as the Normal Forms the- 
orem [28|, we indeed find that the center manifold of E 
attracts an open set of orbits provided A2 < 0. On the 
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Cr. P. 


Xc 


yc 


Zc 


Existence 


Stable manifold 




Wtot 


Acceleration 


D 


Xc4 


yc4 




Always 


l-Dimensional 








Never 


E 


Xc5 




Zc5 


Always 


2-Diniensional 


1 


-1 


Always 



TABLE II: The real and physically meaningful critical points of Model 2 and their behavior. 



other hand, if A2 > the orbits located near the center 
manifold of E blow up in a finite time. Since this point 
does not allow for a solution of the coincidence problem 
(it always possesses ft 4, — 1) we do not present the afore- 
mentioned procedure in detail. 

Numerical investigation reveals the above features. In 
fig. [3] we depict orbits projected in the xy-plane, as they 
arise from numerical evolution in the case of A2 = —0.5 
and /i2 = 0.5. 




C. Model 3: Power-law potential and 
exponentially-dependent dark-matter particle mass 



In this case the autonomous system reads: 



x' = -3x + ^x{l -x^ - y^) - 



,3 2 2n ^■ixyz 

y = i^yy^-x -y ) — 

z' = -xz^. (32) 
The real and physically meaningful critical points are 



(a;c6 = 0, ycfs = 1, 2c6 = 0) , 



Ml 



3 

27r 



Vcl = - TT^, ZcT = 1 , (33) 



FIG. 3: (Color Online) sy-projection 0/ the phase-space of 
Model 2, for the parameter values A2 — —0.5 and fi2 = 0.5. 
The critical point E (representing de Sitter solutions) is the 
attractor of the system. The dashed (red) curves bound the 



and the necessary conditions for their existence are shown 



physical part of the phase space, that ts corresponding to in table UTTl The 3 X 3 matrix Q of the linearized pertur- 



< < 1 



bation equations writes: 



^ j {-9x1 - 2V6^ilXc - 3 (2/2 + 1)) 

Q = -^yd^Xc + ZcXt) 



yc (-3a;c - Zc\2 + \/6^i) 

(-92/2 - Xc{iXa + ZcA2) + 3) 





2 

hxcycM 



In the model at hand, all critical points are non- 
hyperbolic and the dimensionality of their stable 
manifold is presented in table Hill Additionally, we men- 
tion that there exists also an unphysical critical point 

H, with coordinates ixcs, = ~\/ f/^ii 2/c8 = Zcs = Gi- 



lts stable manifold is 2D if > yl, and ID if 

< For the choice \fii\ > the orbits initially 
below the stable manifold of G converge towards H. 
The orbits initially above this curve converge towards 
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Xc 


yc 


Zc 


Existence 


Stable manifold 




Wtot 


Acceleration 


F 


Xc6 


yc6 




Always 


2-Dimensional 


1 


-1 


Always 


G 


Xc7 


yc7 


Zc7 


Iml > ^ 


1-Dimensional 


1-4 


-1 


Always 



TABLE III: The real and physically meaningful critical points of Model 3 and their behavior. 



F. This behavior is depicted in fig. 31 which has arisen 
from numerical evolution using A2 = 1 and fii — 1.8. 



If we restrict ourselves in the region < then 




FIG. 4: (Color Onlme)xy-projecUon of the phase-space of 
Model 3 for the parameter values A2 = 1 and /ii — 1.8. The 
stable manifold of G (thick curve) divides the physical part of 
the phase space (region bounded by the dashed (red) curves) 
in two regions. The orbits initially below this curve converge 
towards H, while those initially above this curve converge to- 
wards F. 



the critical point G does not exists and thus there 



are not scaling solutions. In this case F is indeed the 
attractor for a positive-measure set of initial conditions. 
Moreover, there exist exceptional orbits contained on 
a ID center manifold of F whose dynamical behavior 
cannot be anticipated from the linear analysis. However, 
since this scenario does not lead to a solution of the 
coincidence problem {Q^ = 1 always) we do not present 
an advanced stability analysis for F. 

D. Model 4: Exponential potential and 
power-law-dependent dark-matter particle mass 

In this case the autonomous system writes: 



x' = -Sx + ^x{l - - y2) _ y 

Ai2 



— ^z(l + .T -y ) 



y' = -x^ - 2/^) - Y xy 



z = —xz 



(34) 



The real and physically meaningful critical points are 

{Xc9 = 0, yc9 = 0, Zc9 = 0) , 

(xcw = 2^cio = y^l + y, = j , (35) 



and in table IIVI we present the necessary conditions for 
their existence. The 3x3 matrix Q of the linearized 
perturbation equations reads: 



i (-9a;2 _ 2z^ii2Xc ~ "iyl ~ 3) 

Q= -i2/c(6xe + V6A1) 

.,2 



yc {-ixc - \/6Ai + Zcti2) 
-9yl - Xc (3xe + VeXi) + 3) 








1 M2 



The aforementioned critical points are non-hyperbolic 
since at least one eigenvalue of Q is always zero. Lin- 
ear analysis in not conclusive in these cases, but infor- 
mation about the dimensionality of the stable manifold 
can be obtained by applying the center manifold theorem 
[2^ . The corresponding results are shown in table IIVI 
Since both I and J cannot solve the coincidence problem 



(510 — 1), we do not present the aforementioned analysis 
in detail. Finally, in order to acquire a more transparent 
picture of the phase-space behavior, we evolve the system 
numerically for Ai = 1 and fj,2 = 1-8 and we depict the 
results in fig. [51 
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Xc 


yc 


Zc 


Existence 


Stable manifold 




Wtot 


Acceleration 


I 


Xc9 




Zc9 


Always 


1-Dimensional 








Never 


J 






ZclO 


Always 


2-Dimensional 


1 


-i(3 + A?) 


Always 



TABLE IV: The real and physically meaningful critical points of Model 4 and their behavior. 
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FIG. 5: (Color Online) i:?/-projecfion of the phase-space of 
Model 4 for the parameter values Ai = 1 and ^2 = 1.8. 
The critical point J (corresponding to a super- accelerating um- 
verse) attracts all the orbits in this invariant set. The dashed 
(red) curves bound the physical part of the phase space, that 
is corresponding to < fi^ < 1. 



IV. COSMOLOGICAL IMPLICATIONS AND 
DISCUSSION 

Having performed a complete phase-space analysis of 
several varying dark-matter-mass models, we can discuss 
the corresponding cosmological behavior. A general re- 
mark is that this behavior is radically different from the 
corresponding quintessence scenarios with the same po- 
tentials and mass-functions 0, [li, [H, [13, d, El, US [IH . 
Additionally, a common feature of almost all the phan- 
tom models previously studied is the existence of attrac- 
tors with W0 < — 1 in the whole phase-space and 
thus, independently of the specific scenario and of the 
imposed initial conditions, the universe always lies below 
the phantom divide, as it is expected for phantom cos- 
mology. This global behavior is not always realized in 
the case of exponentially dependent dark-matter mass, 
and additional constraints must be imposed. 

Apart form acquiring acceleration, in this work we 
examine whether the above constructed varying dark- 
matter-mass models can solve or alleviate the coincidence 
problem. Thus, assuming as usual that the present uni- 
verse is already at a late-time attractor, we calculate 



in all stable fixed points, and if < Sl^ < 1 then the 
coincidence problem is solved since fi^ and ^dm will be 
of the same order of magnitude as suggested by observa- 
tions. On the contrary, = 1 corresponds to a universe 
completely dominated by dark energy, while fJ^ = (that 
is n DM = 1 ) to one completely dominated by dark mat- 
ter, both in contrast with observations. 

Finally, we mention that as long as the interaction re- 
sponsible for the varying dark-matter particle mass is 
not too strong, the standard cosmology can be always 
recovered. On the other hand, since we assume that the 
universe is currently at an attractor, its state is indepen- 
dent of the initial conditions. Thus, we can switch on 
the interaction and consider as initial conditions the end 
of the known epochs of standard Big Bang cosmology, in 
order to avoid disastrous interference. 



A. Model 1 

In this model the critical point B is unstable, and there- 
fore it cannot be a late-time cosmological solution. The 
only relevant critical point is A, which is a stable fixed 
point for Ai (/ii — Ai) < 3. As can be seen from tablelH it 
corresponds to an accelerating universe with fl^ = 1, that 
is to complete dark-energy domination. Thus, this spe- 
cific cosmological solution cannot solve the coincidence 
problem. Furthermore, the fact that Wtot is not only 
less than —1/3, as required by the acceleration condi- 
tion, but it is always less than —1, leads to ij > at all 
times. Therefore, this solution corresponds to a super- 
accelerating universe (30j . that is with a permanently 
increasing H, resulting to a Big Rip. This behavior is 
common in phantom cosmology [j, |3l| . 

A remarkable feature of this model, as well as of Model 
3, is that if there exist scaling solutions, then, for a wide 
region of the parameter space, the stable manifold of the 
corresponding critical point marks the basin of attraction 
of either a phantom attractor or an unphysical attracting 
state. Thus, there exist an open set of orbits of the physi- 
cal region that converge to an unphysical state instead to 
a phantom solution. This behavior was revealed analyt- 
ically and it was confirmed through numerical elabora- 
tion, and seems to be typical in the case of exponentially- 
dependent dark-matter mass in the phantom framework. 
To avoid dealing with unphysical states, we can either 
restrict the physical portion of the phase-space to the re- 
gion above the stable manifold of the scaling solutions, 
or restrict the parameter-space itself. In both cases we 
obtain an additional constraint, that was not present in 
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previous studies of phantom cosmology [1, [T^ , which 
further weakens the apphcabihty of the model. 

In summary, Model 1, that is an exponential poten- 
tial and an exponentially-dependent dark-matter particle 
mass, cannot act as a candidate for solving the coinci- 
dence problem. 



B. Model 2 

In this case, both real and physically meaningful crit- 
ical points, namely D and E, have a stable manifold of 
smaller dimensionality than that of the phase-space, and 
as was mentioned in subsection llll BI thev have very small 
probability to be the late-time attractors of the system. 
However, even if the cosmological evolution is managed 
to be attracted by these solutions, the coincidence prob- 
lem will not be solved, since D represents a flat, non- 
accelerating universe dominated by dark matter, and E 
correspond to de Sitter universe completely dominated 
by dark energy. These critical points are located in the 
region where the scalar field and the Hubble parameter 
diverge. Divergencies in a cosmological scenario are rep- 
resented as asymptotic states, in particular associated 
with the past and future asymptotic dynamics pU . Is^ l . 
In the present Model 2, due to the non-compactness of 
the phase-space, such a behavior can lead either to an 
asymptotic state acquired at infinite time, or to a singu- 
larity reached at a finite time. If if — > oo or — + oo at 
t oo then we acquire an eternally expanding universe, 
while if iJ — !■ cxD at t ^ Ibr < oo then the universe 
results to a Big Rip [33l |. 

Therefore, power-law potentials with power-law- 
dependent dark-matter particle masses, cannot solve the 
coincidence problem. 



C. Model 3 

In this model we see that the critical point F exists 
always, while G exists only for \fii\ > V3- However, in 
both cases the stable manifold is of smaller dimensional- 
ity than that of the phase-space. Furthermore, in order 
to avoid the treatment of unphysical attracting states 

we have to impose the additional constraint < 

For this choice of parameters, G does not exists and thus 
there are not scaling solutions, while F is the attractor for 
a positive-measure set of initial conditions. Point F cor- 
responds to a dark-energy dominated de Sitter universe, 
while G to a flat accelerating universe with — 1 — , 

that is with < fi^ < 1 in the region that it exists. In 
both points the phantom field diverges. However, even if 
G possesses < Sl^ < 1, it can not solve the coincidence 
problem since it is not a relevant late-time attractor. 

In summary, power-law potentials with exponentially- 
dependent dark-matter particle masses cannot solve or 
even alleviate the coincidence problem. 



D. Model 4 

In this case, the critical points I and J exist al- 
ways. The point I corresponds to a flat, non-accelerating, 
matter-dominated universe. J corresponds to a dark- 
energy dominated universe, that super-accelerates [30[. 
However, similarly to the previous cases, the stable mani- 
folds of I and J are ID or 2D respectively, and thus almost 
all orbits of the cosmological system cannot be attracted 
by them at late times. In addition, they cannot lead to 
< < 1. Therefore, an exponential potential and a 
power-law-dependent dark-matter particle mass, cannot 
solve the coincidence problem. 



V. CONCLUSIONS 

In the present work we investigated the phantom cos- 
mological scenario, with varying-mass dark-matter parti- 
cles due to the interaction between dark-matter and dark- 
energy sectors. In particular, we performed a detailed 
phase-space analysis of various models, with either ex- 
ponentially or power-law dependent dark-matter particle 
mass, in exponential or power-law scalar field potentials. 
These functions cover a wide range of the possible forms, 
and they correspond to the cases that can accept a rea- 
sonable theoretical justification [H, [H, [11, [13, El, [l^l . In 
each case we extracted the critical points, we determined 
their stability, and we calculated the basic cosmologi- 
cal observables, namely the total equation-of-state pa- 
rameter wtot and r^DE (attributed to the phantom field). 
Our basic goal was to examine whether there exist late- 
time attractors, corresponding to accelerating universe 
and possessing ^be/^dm ~ C'Cl); thus satisfying the 
basic observational requirements. 

Apart from the case of an exponential potential with 
an exponentially-dependent dark-matter particle mass, 
which possesses a relevant late-time (phantom) attrac- 
tor, in all the other models we found that physical, well- 
motivated solutions have a very small chance to attract 
the universe at late times. In addition, in all the exam- 
ined cases, solutions having JIde/^^dm ~ C'(l) are not 
relevant attractors at late times. Therefore, summariz- 
ing, the coincidence problem cannot be solved or even 
alleviated in varying-mass dark matter particles mod- 
els in the framework of phantom cosmology, in a rad- 
ical contrast with the corresponding quintessence case 
[H, S HI- This conclusion agrees with that of [1], 
that interacting phantom cosmology cannot solve the co- 
incidence problem. It seems that interacting phantom 
cosmology, either directly or through the dependence of 
the dark-matter particle mass, cannot fulfill the basic re- 
quirements that led to its construction, that is to provide 
stable accelerating late-time solutions which can solve the 
coincidence problem. An alternative direction could be to 
consider a specially constructed potential or dark-matter 
particle mass in order to solve the coincidence problem, 
but this would imply significant loss of simplicity, gener- 
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ality, and theoretical justification of the model. 

The aforementioned conclusion has been extracted 
by the negative-kinetic-cnergy realization of phantom, 
which does not cover the whole class of phantom mod- 
els. However, since it is a qualitative statement it should 
intuitively be robust for general phantom scenarios, too. 
Therefore, phantom cosmology with varying-mass dark 
matter particles cannot easily act as a successful candi- 



date to describe dark energy. 
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